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Neutrino telescopes that measure relative fluxes of ultrahigh-energy v e ,Vn, v T can give information 
about the location and characteristics of sources, about neutrino mixing, and can test for neutrino 
instability and for departures from CPT invariance in the neutrino sector. We investigate conse- 
quences of neutrino mixing for the neutrino flux arriving at Earth, and consider how terrestrial 
measurements can characterize distant sources. We contrast mixtures that arise from neutrino os- 
cillations with those signaling neutrino decays. We stress the importance of measuring v e ,v )i ,v T 
. fluxes in neutrino observatories. 

o 
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Ch ! I. INTRODUCTION 

Neutrino telescopes promise to probe the deepest reaches of stars, galaxies, and exotic structures in the cosmos 0, 
£S) | H,H,I3' Unlike charged particles, neutrinos arrive on a direct line from their source, undefiected by magnetic fields. 

Unlike photons, neutrinos interact weakly, so they can penetrate thick columns of matter 20]. It is plausible that 
ultrahigh-energy extraterrestrial neutrinos will emerge from the atmospheric-neutrino background at energies between 
1 and 10 TeV. Prospecting for extraterrestrial neutrino sources leads the science agenda for neutrino observatories, 
with characterizing the sources and the processes that operate within them to follow. The fluxes of extraterrestrial 
neutrinos may, in addition, offer clues to the properties of neutrinos themselves. 

The task of developing diagnostics for neutrino sources by measuring relative fluxes of electron-, muon-, and tau- 
neutrinos is complicated by neutrino oscillations. Prominent among expected sources is the diffuse flux of neutrinos 
produced in the jets of active galactic nuclei such as the TeV-gamma sources Mkn 421 and 501, which are some 
140 Mpc distant from Earth. The vacuum oscillation length, L osc — 4ffE v /\Am 2 \, is short compared with such 
intergalactic distances, so neutrinos oscillate many times between source and detector. For |Am 2 | = 10 -5 eV 2 , for 
example, the oscillation length is L osc sa 2.5 x 10~ 24 Mpc • {E v /1 eV), a fraction of a megaparsec even for 10 20 -eV 
neutrinos. The fluxes $ = {ip e , (p^, (p T } that arrive at Earth are not identical to the source fluxes <I> = {tp% y° , 
and the transfer matrix X that maps <I> to $ is not, in general, invertible |2l|. Moreover, over their long flight 
paths, cosmic neutrinos are vulnerable to decay processes that would have gone undetected in terrestrial or solar 
experiments. Should the neutrino sector exhibit departures from CPT invariance, the pattern of neutrino mixing 
could be considerably richer — and less circumscribed by experiment — than conventional wisdom holds. 

In this paper, we investigate the consequences of neutrino mixing for the cosmic neutrino flux $ at Earth, and explore 
how measurements on Earth can characterize the source flux We contrast the fluxes that result from neutrino 
mixing with those that might arise from simple decay scenarios, and we look at the unconventional consequences that 
might obtain if CPT symmetry were violated. We first carry out an idealized analysis, taking our cue from current 
experiments; then we take into account the uncertainties of existing experimental constraints; finally we ask what 
we will know after the next round of neutrino oscillation experiments. Our aim is to examine the scientific potential 
of neutrino observatories in light of current information about neutrino properties, and to project the situation five 
years hence. 

The basic operational goal of neutrino telescopes is the detection of energetic — hence, long-range — muons generated 
in charged-current interactions {v^D^N — > + anything. Efficient, well-calibrated detection of {v e ,v e ) 

interactions is also required to make neutrino observatories incisive tools for the investigation of cosmic sources 
and neutrino properties. Good detection of (y T ,v T ) interactions would test the expectation that muon neutrinos 
and tau neutrinos arrive at Earth in nearly equal numbers. The normalized electron-neutrino flux at Earth emerges 
as a promising diagnostic for the character of cosmic sources and for nonstandard neutrino properties. If neutrinos 
behave as expected — mixing according to the standard three-generation picture — it should prove possible, over time, 
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to infer something about the flavor mix of neutrinos produced in distant sources. As is by now well known, neutrino 
oscillations — conventionally understood — map the standard neutrino mixture at the source, <£>j? td = {A,|,0} into 
a mixture at Earth that approximates {|, |, | }. We explore the uncertainties that attach to this expectation in 
light of recent improvements in the experimental constraints on neutrino- mixing parameters, and we compute the 
mixture at Earth to be expected for other source mixtures. Neutrino decays over the long path from astrophysical 
sources could distort the flavor mixture of neutrinos arriving at Earth. Current speculations about CPT violation in 
the neutrino sector imply small, and perhaps undetectable, modifications to the conventional oscillation scenario for 
ultrahigh-energy neutrinos, but would lead to striking consequences for antineutrino decay. 



II. THE INFLUENCE OF NEUTRINO OSCILLATIONS 



Neutrino emission from active galactic nuclei (AGNs) may constitute the dominant diffuse flux at energies above a 
few TeV, where cosmic sources should emerge from the background of atmospheric neutrinos. With luminosities on 
the order of I0 45±3 erg/s, AGNs are the most powerful radiation sources in the universe. They are cosmic accelerators 
powered by the gravitational energy of matter falling in upon a supermassive black hole. Protons accelerated to very 
high energies within an AGN may interact with ultraviolet photons in the bright jets along the rotation axis or with 
matter in the accretion disk. The resulting pp or pj collisions yield approximately equal numbers of 7r + , tt°, tt~ that 
decay into /i + z^, 77, or pTv^. The subsequent muon decays yield additional muon neutrinos and electron neutrinos, 
so the final count from 7r + + 7r° + 7r~ is 27 + 2v^ + 2v^ + lv e + lD e , for a normalized neutrino flux 

^ = {^ = 1,^ = 1,^=0} (1) 

at the source. This standard mixture — equally divided between neutrinos and antineutrinos — is the canonical expec- 
tation from most hypothesized sources of ultrahigh-energy neutrinos. 

In the standard, three-generation description of neutrino oscillations, the flavor eigenstates v a are related to the 
mass eigenstates Vi through v a = ^ UpiVi, or 
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A canonical form for the neutrino mixing matrix is p| 



/ U e i U e2 U e3 \ / C12C13 S12C13 si 3 e~ lS \ 

U = I Uf+i 17^2 Uf+3 J = I -S12C23 - Ci2S23Sl3e 1 ' 5 C12C23 - Sl2S23Sl3e 115 S23C13 J ) (3) 

\U T i U t2 U t3 J \ s 12 s 23 - c 12 c 2 3Si 3 e lS -c 12 s 23 - s 12 c 23 Si 3 e lS c 23 c 13 j 

where Sij = sin^j, cy = cos#y and S is a CP-violating phase. A useful idealized form of the mixing matrix that 
holds for sin 0i3 = 0, and maximal atmospheric mixing, sin 2#23 = 1, is 

(C12 S12 \ 

S12/V2 C12/V2 1/V2 . (4) 
S12/V2 -C12/V2 1/V2J 

We can express the oscillation (or survival) probabilities in a transfer matrix X that maps an initial mixture of 
flavors at the source into the observed mixture at the detector, 




(5) 

Averaged over many oscillations, X takes the form 

Xpa = S a — 2jfty^ UgjUpiUgjUg-j 

i>j 

= E 1 u * 1 2 1 u n 1 2 . ( 6 ) 
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FIG. 1: Ternary plots of the neutrino flux $ at Earth, showing the implications of current (left pane) and future (right pane) 
knowledge of neutrino mixing. The small black squares indicate the v R fractions produced by the idealized transfer matrix 
"^ideai as ip^ varies from to 1 in steps of 0.1. A crossed circle marks the standard mixed spectrum at Earth, <E> st d = {|, g, ^}; 
for comparison, a red dot marks the standard source spectrum, $g td = {^, |,0}. Colored swaths delimit the fluxes at Earth 
produced by neutrino oscillations from the source mixtures <E>g = {0,1,0} (pink), $° td (red), and = {1,0,0} (orange), 
using 95% CL ranges for the oscillation parameters. Black crosses (x) show the mixtures at Earth that follow from neutrino 
decay, assuming normal (ip e ~ 0.7) and inverted (<p e ~ 0) mass hierarchies. The blue bands show how current and future 
uncertainties blur the predictions for neutrino decays. The violet tripod indicates how CPT-violating oscillations shape the mix 
of antineutrinos that originate in a standard source mixture, and the violet cross averages that i> mixture with the standard 
neutrino mixture. The brown squares denote consequences of CPT violation for antineutrino decays. 



which does not depend on the phase 5. The idealized form that corresponds to is 

(1 — 2x x x \ 

x Ul-x) \{l-x) , (7) 
x |(1 -x) |(1 - x) ) 

where x = sin 2 12 cos 2 12 . Because the second and third rows are identical, the and v T fluxes that result from 
any source mixture $° are equal: ip^ — tp T . 

It is noteworthy that Aidoai maps the standard source mixture $!? td into 

$std = {Ve = |,V/t = f^T = §}) (8) 

independent of x, i.e., independent of #i 2 . Because and v T are fully mixed, the result is more general: any source 
with <p® = i is mapped into $ s td, independent of x. The source mixture = {5, §, §} that arises from Z° decay is yet 
more special because it is unchanged by neutrino oscillations: any three- generation transfer matrix maps it into ^std- 
For an arbitrary value of the v e fraction at the source, X- X &<>a\ leads to <p e = — 3x) + x, with <p M = ip T = |(1 — <p e ). 
The variation of (p e with the v e source fraction <p® is shown as a sequence of small black squares (for tp° = 0, 0.1, . . . ,1) 
in FigureHfor the value x — 0.21, which corresponds to 12 = 0.57, the central value in a recent global analysis [||. 
The v e fraction at Earth ranges from 0.21, for (p® — 0, to 0.59, for ip® = 1. 

The simple analysis based on Aideai is useful for orientation, but it is important to explore the range of expectations 
implied by global fits to neutrino-mixing parameters. The recent KamLAND data 7] on the reactor-j/ e rate and 
spectral shape, taken together with solar-neutrino data, select the large-mixing- angle (LMA) solution, with a mixing 
angle bounded at 95% CL to lie in the interval 0.49 < 612 < 0.67 6]. Observations of atmospheric neutrinos favor 
maximal mixing; we take j x 0.8 < 823 < ? X 1.2 @, also at 95% CL. Experimental evidence favors #13 <C 1; informed 
by the nonobservation of oscillations in the CHOOZ and Palo Verde reactor experiments |9lll0|. we take < 6*13 < 0.1. 

With current uncertainties in the oscillation parameters, a standard source spectrum, <J>!? td = {^, |,0}, is mapped 
by oscillations onto the red boomerang near <J> st d = {5, 3} in the left pane of Figure ^ Given that A^cai maps 
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$g td — > $ st d for any value of 0i2, it does not come as a great surprise that the target region is of limited extent, and 
this is a familiar result [1 lj . The variation of #23 away from ^ breaks the identity ip^ = (p T of the idealized analysis. 

One goal of neutrino observatories will be to characterize cosmic sources by determining the source mix of neutrino 
flavors. It is therefore of interest to examine the outcome of different assumptions about the source. We show in 
the left pane of Figure ^ the mixtures at Earth implied by current knowledge of the oscillation parameters for source 
fluxes $[] = {0, 1, 0} (the purple band near tp e ps 0.2) and $5 = {1) 0, 0} (the orange band near (p e ps 0.6) j^. For the 
$g td and $5 source spectra, the uncertainty in 9\2 is reflected mainly in the variation of ip e , whereas the uncertainty 
in #23 is expressed in the variation of tp^/tpr For the $g case, the influence of the two angles is not so orthogonal. 
For all the source spectra we consider, the uncertainty in 813 has little effect on the flux at Earth. The extent of the 
three regions, and the absence of a clean separation between the regions reached from <&° td and $g indicates that 
characterizing the source flux will be challenging, in view of the current uncertainties of the oscillation parameters. 

Over the next five years — roughly the time scale on which large-volume neutrino telescopes will come into 
operation — we can anticipate improved information on 9jo and #23 from KamLAND |l3j| and the long-baseline ac- 
celerator experiments at Soudan ^ij and Gran Sasso We base our projections for the future on the ranges 
0.54 < 6>i2 < 0.63 and J x 0.9 < #23 < f x 1.1, still with < 9 13 < 0.1 The results are shown in the right panel of 
Figure n The (purple) target region for the source flux $q shrinks appreciably and separates from the (red) region 
populated by $g td , which is now tightly confined around ^std- The (orange) region mapped from the source flux $5 
by oscillations shrinks by about a factor of two in the ip e and ip„, — ip T dimensions. 

If CPT invariance is not respected in the neutrino sector |l6j, the neutrino and antineutrino mixing matrices need 
not coincide. The KamLAND experiment's confirmation of solar-like oscillations of v e restricts the space of CPT- 
violating mixing matrices, but Barenboim, Borissov, and Lykken [l7j have found two illustrative solutions that respect 
all existing constraints: 

/ 012 = 0.6 \ /0i2 = 0.785 \ 

Solution 1: 23 = 0.5 ; Solution 2: 6» 23 = 0.52 . (9) 
\0i3 = O.O8/ Wi 3 = 0.08/ 

These two solutions lead to essentially identical transfer matrices (cf. Eqn. ©) that map a standard source mixture of 
animeutrinos into $cptv = {0.42, 0.36, 0.22}, which is plotted as a violet tripod in Figure^ Observing <p e /(p e = § 7^ 1 
would be manifest evidence of CPT violation. Although in principle it might be possible to distinguish the flux of v e 
and v e at Earth by observing the resonant formation v e e — > W~ — > hadrons, the v e + v e flux is more likely to be the 
available observable. In the presence of CPT-violating mixing, the expectation from a standard flux at the source is 
($Cptv) = |(^std + ^cptv) = {0.375,0.35,0.275}, which is plotted as a violet cross in Figure^ Distinguishing this 
value from $ st d = {^, \ 1 5} would be extraordinarily challenging. 



III. RECONSTRUCTING THE NEUTRINO MIXTURE AT THE SOURCE 



What can observations of the blend $ of neutrinos arriving at Earth tell us about the source? Inferring the nature 
of the processes that generate cosmic neutrinos is more complicated than it would be in the absence of neutrino 
oscillations. Because and v r are fully mixed — and thus enter identically in Aid ea i — it is not possible fully to 
characterize $°. We can, however, reconstruct the v e fraction at the source as 

The reconstructed source flux <p° e is shown in Figure [21 as a function of the v e flux at Earth. The heavy solid line 
represents the best- fit value for 0i2; the light blue lines and thin solid lines indicate the current and future 95% CL 
bounds on i2 . 

A possible strategy for beginning to characterize a source of cosmic neutrinos might proceed by measuring the v^Jv^ 
ratio and estimating ip e under the plausible assumption — later to be checked — that ip^ = ip T . Let us first note that 
very large (cp e > 0.65) or very small (cp e < 0.15) v e fluxes cannot be accommodated in the standard neutrino-oscillation 
picture. Observation of an extreme v e fraction would implicate unconventional physics. 

As we have already seen from Figure^ constraining the source flux sufficiently to test the nature of the neutrino 
production process will require rather precise determinations of the neutrino flux at Earth. Suppose we want to test 
the idea that the source flux has the standard composition of Eqn. QJ. With today's uncertainty on 12 , a 30% 
measurement that locates <p e = 0.33 ± 0.10 implies only that < (p® < 0.68. For a measured flux in the neighborhood 
of i, the uncertainty in the solar mixing angle is of little consequence: the constraint that arises if we assume the 
central value of 0i2 is not markedly better: 0.06 < pP e < 0.59. A 10% measurement of the v e fraction, <p e — 0.33±0.033, 



0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

<l>e 



FIG. 2: The source flux 95° of electron neutrinos reconstructed from the u e flux <p e at Earth, using the ideal transfer matrix 
^ideai of Eqn. The heavy solid line refers to our chosen central value, #12 = 0.57. The light blue lines refer to the current 
experimental constraints (at 95% CL), and the thin solid lines refer to our projection of future experimental constraints. 

would make possible a rather restrictive constraint on the nature of the source. The central value for 6*12 leads to 
0.26 < ifP e < 0.43, blurred to 0.22 < < 0.45 with current uncertainties. 

Measured v e fractions that depart significantly from the canonical ip e = ^ would suggest nonstandard neutrino 



sources. An observed flux ip e — 0.5 ± 0.1 points to a source flux 0.47 < cp® < 1, with current uncertainties, whereas 
<p e = 0.25 ± 0.10 indicates < ip° e < 0.35. 



Beacom, Bell, and collaborators LL9j have observed that the decays of unstable neutrinos over cosmic distances 
can lead to mixtures at Earth that are incompatible with the oscillations of stable neutrinos. The candidate decays 
are transitions between mass eigenstates by emission of a very light particle, 2/j — * {v rj ,Dj) + X. Dramatic effects 
occur when the decaying neutrinos disappear, either by decay to invisible products or by decay into active neutrino 
species so degraded in energy that they contribute negligibly to the total flux at the lower energy. If the lifetimes of 
the unstable mass eigenstates are short compared with the flight time from source to Earth, decay of the unstable 
neutrinos will be complete, and the (unnormalizcd) flavor v a flux at Earth is given by 



with ip a = (fia/Y^pVp- 

If only the lightest neutrino survives, the flavor mix of neutrinos arriving at Earth is determined by the flavor 
composition of the lightest mass eigenstate, independent of the flavor mix at the source. (The flavor mix at the source 
of course determines the relative fluxes of the mass eigenstates at the source, and so the rate of surviving neutrinos at 
Earth.) For a normal mass hierarchy mi < 1712 < m^, the v a flux at Earth is </? Q = |J7 a i| 2 . Accordingly, the neutrino 
flux at Earth is $ norma ] = {|J7 e i| 2 , |£^»i| 2 , |^ri| 2 } ~ {0.70, 0.17, 0.13} for our chosen central values of the mixing angles. 
If the mass hierarchy is inverted, mi > m2 > TO3, the lightest (hence, stable) neutrino is ^3, so the flavor mix at Earth is 
determined by ip a — |J7 Q 3| 2 . In this case, the neutrino flux at Earth is ^inverted = {|^e3| 2 , |^3| 2 > | | 2 } ~ {0, 0.5, 0.5}. 
Both ^normal and ^inverted, which are indicated by crosses (x) in Figure ^ are very different from the standard flux 
<I> s td = {fe = 5, <Ptt = = \} produced by the ideal transfer matrix from a standard source. Observing either 

mixture would represent a departure from conventional expectations. 

The fluxes that result from neutrino decays en route from the sources to Earth are subject to uncertainties in the 
neutrino- mixing matrix. The expectations for the two decay scenarios are indicated by the blue regions in Figure ^ 



IV. INFLUENCE OF NEUTRINO DECAYS 
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where we indicate the consequences of 95% CL ranges of the mixing parameters now and in the future. With current 
uncertainties, the normal hierarchy populates 0.61 < ip e < 0.77, and allows considerable departures from tp^ = ip T . 
The normal-hierarchy decay region based on current knowledge overlaps the flavor mixtures that oscillations produce 
in a pure-z/ e source, shown in orange. (It is, however, far removed from the standard region that encompasses 
^std-) With the projected tighter constraints on the mixing angles, the range in ip e swept out by oscillation from a 
pure-Ve source or decay from a normal hierarchy shrinks by about a factor of two. Neutrino decay then populates 
0.65 < (fi e < 0.74, and is separated from the oscillations. The degree of separation between the region populated by 
normal-hierarchy decay and the one populated by mixing from a pure-^ e source depends on the value of the solar 
mixing angle 9i2- For the seemingly unlikely value #12 = § , both mechanisms yield $ = {|, |, |}. 

With current uncertainties, the inverted hierarchy spans the range {(p^.Lpr) — (0.34,0.66) to (0.66,0.34), always 
with ip e ps 0. With the improvements we project in the knowledge of mixing angles, the range will decrease to 
(ip^jipr) = (0.42,0.58) to (0.58,0.42), with <p e w 0. In both cases, these mixtures are well separated from the 
mixtures that would result from neutrino oscillations, for any conceivable source at cosmic distances. 

Should CPT not be a good symmetry for neutrinos, the mass hierarchies and mixing matrices can be different 
for neutrinos and antineutrinos. Different numbers of stable v and v may reach the Earth, so we do not know how 
to combine neutrino and antineutrino results to compare with neutrino-telescope measurements. The analysis of 
neutrino decays does not change if we relax CPT invariance, because neutrino mixing is already highly constrained by 
experiment. The Barenboim-Borissov-Lykkcn antineutrino mixing parameters 17] (cf. Eqn. @) provide examples 
of what can be expected for antineutrino decays if CPT is violated. Consider first the normal mass hierarchy, 
in which v-i and V3 both decay completely before reaching Earth. The resulting antineutrino fluxes at Earth are 
§Cptvi = {0.68,0.28,0.04} and <E>cptv2 = {0.50,0.42,0.08}, both characterized by unusually large ratios of (p^/pr 
that place them outside the range of CPT-conserving decays and oscillations alike j2^|. In the case of an inverted 
mass hierarchy, both BBL parameter sets predict $cptv inv = {0.01,0.23,0.76}, an unusually small ratio of tp^/ipr. 
Conventional oscillations — or CPT-conserving decays — populate the ternary plots of Figure ^ along the line ps^ — tp T . 
CPT-violating decays offer an example of exotic physics with a very specific signature. 

Beacom and Bell [l^ have shown that observations of solar neutrinos set the most stringent plausible lower bound 
on the reduced lifetime of a neutrino of mass m as r/m > 10~ 4 s/cV. This rather modest limit opens the possibility 
that some neutrinos do not survive the journey from astrophysical sources, with consequences we have just explored. 
The energies of neutrinos that may be detected in the future from AGNs and other cosmic sources range over several 
orders of magnitude, whereas the distances to such sources vary over perhaps one order of magnitude. The neutrino 
energy sets the neutrino lifetime in the laboratory frame; more energetic neutrinos survive over longer flight paths than 
their lower-energy companions 24]. Under propitious circumstances of reduced lifetime, path length, and neutrino 
energy it might be possible to observe the transition from more energetic survivor neutrinos to less energetic decayed 
neutrinos. 

If decay is not complete, the (unnormalized) flavor v a flux arriving at Earth from a source at distance L is given by 



with the normalized flux p a (E v ) = p a (E v )/ <pp(E v ). An idealized case will illustrate the possibilities for observing 
the onset of neutrino decay and estimating the reduced lifetime. Assume a normal mass hierarchy, mi < m,2 < ITI3, 
and let T3/TO3 = T2/TO2 = r/m. For a given path length L, the neutrino energy at which the transition occurs 
from negligible decays to complete decays is determined by r/m. We show in the left pane of Figure |21 the energy 
evolution of the normalized neutrino fluxes arriving from a standard source; the energy scale is appropriate for 
the case r/m = 1 s/eV and L = 100 Mpc. The result shown there is actually universal — within our simplifying 
assumptions — provided that the neutrino energy scale is renormalized as 



We show in the left pane of Figure 0] that over about one decade in energy, the flavor mix of neutrinos arriving at 
Earth changes from the normal decay flux, ^normal, to the standard oscillation flux, <I> st d. For the case at hand, the 
transition takes place between 3 x 10 15 eV and 3 x 10 16 eV. The corresponding results for the inverted mass hierarchy 
are depicted in the right panes of Figures and ^ which show a rapid transition from "^inverted to "I'std- 

If we locate the transition from survivors to decays at neutrino energy E* , then we can estimate the reduced lifetime 
in terms of the distance to the source as 




(12) 




(13) 




(14) 
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FIG. 3: Energy dependence of normalized v e , Vy,, and v T fluxes, for the two-body decay of the two upper mass eigenstates, 
with the neutrino source at L = 100 Mpc from Earth and r/m = 1 s/eV. The left pane shows the result for a normal mass 
hierarchy; the right pane shows the result for an inverted mass hierarchy. With suitable rescaling of the neutrino energy (cf. 
Eqn. 11310 . these plots apply for any combination of path length and reduced lifetime. 

In practice, ultrahigh-energy neutrinos are likel y t o arrive from a multitude of sources at different distances from 
Earth, so the transition region will be blurred |2q . Nevertheless, it would be rewarding to observe the decay-to- 
survival transition, and to use Eqn. I|14|l to estimate — even within one or two orders of magnitude — the reduced 
lifetime. If no evidence appears for a flavor mix characteristic of neutrino decay, then Eqn. i|14|) provides a lower 
bound on the neutrino lifetime. For that purpose, the advantage falls to large values of L/E*, and so to the lowest 




FIG. 4: Ternary plots showing the energy dependence of normalized u e , v^, and v T fluxes, for the two-body decay of the two 
upper mass eigenstates, with the neutrino source at L = 100 Mpc from Earth and r/m = 1 s/eV. The left pane shows the result 
for a normal mass hierarchy; the right pane shows the result for an inverted mass hierarchy. Each unlabeled step multiplies e 
by 10 1//4 . With suitable rescaling of the neutrino energy (cf. Eqn. 1131 1. these plots apply for any combination of path length 
and reduced lifetime. 
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energies at which neutrinos from distant sources can be observed. Observing the standard flux, <j> st d = {3, i, 3}, 
which is incompatible with neutrino decay, would strengthen the current bound on r/m by some seven orders of 
magnitude, for 10-TeV neutrinos from sources at 100 Mpc pjfl . 

V. ASSESSMENT 

The conventional mechanism for ultrahigh-energy neutrino production in astrophysical sources, when modulated 
by the standard picture of neutrino oscillations, leads to approximately equal fluxes of v e , v^, and v T at a terrestrial 
detector. Observing a strong deviation from this expectation will indicate that we must revise our conception of 
neutrino sources — or that neutrinos behave in unexpected ways on the long journey from source to Earth. With 
the more precise understanding of neutrino mixing that we expect will develop over the next five years, it may 
become possible to characterize the neutrino mixture at the source. On the other hand, detecting equal fluxes of the 
three neutrino flavors can strengthen existing bounds on neutrino lifetimes quite dramatically. If neutrinos do decay 
between source and Earth, the flavor mixture at Earth will be incompatible with the consequences of the standard 
source mixture. As constraints improve on the mixing angles 612 and #23, the flavor mix that results from neutrino 
decay will separate decisively from what can be generated by any source mix plus oscillations. In the best imaginable 
case for neutrino decays, the neutrino lifetime might reveal itself in an energy dependence of the flavor mix observed 
at Earth. If neutrinos decay, the fluxes at Earth are very different for normal and inverted mass hierarchies. 

Conventional neutrino oscillations — and conventional decays — lead to nearly equal fluxes of and v T , and so 
populate only a portion of the Lp e -(p^-Lp T ternary plot. A neutrino mixture at Earth that is far from the (p^ — tp T line 
points to unconventional neutrino physics. CPT violation offers one example of novel behavior, but its predictions 
are tested most cleanly by observing differences between neutrinos and antineutrinos, for which the experimental 
prospects are distinctly limited at ultrahigh energies. 

The analysis presented here offers further motivation to develop techniques for identifying interactions of ultrahigh- 
energy v e , and v T — and for measuring the neutrino energies — in cubic-kilometer-scale neutrino observatories. If 
this can be achieved, prospects for learning about the properties of neutrinos and about the nature of cosmic neutrino 
sources will be greatly enhanced. 

Acknowledgments 

We thank Stephen Parke for a helpful observation, and thank John Beacom and Nicole Bell for lively discussions. 
Fermilab is operated by Universities Research Association Inc. under Contract No. DE-AC02-76CH03000 with the 
U.S. Department of Energy. One of us (C.Q.) is grateful for the hospitality of the Kavli Institute for Theoretical 
Physics during the program, Neutrinos: Data, Cosmos, and Planck Scale. This research was supported in part by the 
National Science Foundation under Grant No. PHY99-07949. 



[1] J. G. Learned and K. Mannheim, Annu. Rev. Nucl. Part. Set. 50, 679-749 (2000). 
[2] F. Halzen and D. Hooper, Kept. Prog. Phys. 65, 1025 (2002). 

[3] F. Halzen, "High-Energy Neutrino Astronomy: Science and First Results," [arXiv:astro-ph/0301143]. 
[4] R. Gandhi, C. Quigg, M. H. Reno, and I. Sarcevic, Astropart. Phys. 5, 81-110 (1996); Phys. Rev. D58, 093009 (1998). 
[5] Peter Fisher, Boris Kayser, and Kevin S. McFarland, Annu. Rev. Nucl. Part. Set. 49, 481-527 (1999). 
[6] P. C. de Holanda and A. Yu. Smirnov, "LMA MSW Solution of the Solar Neutrino Problem and First KamLAND Results," 
[arXiv:hep-ph/0212270]. 

[7] K. Eguchi, et al. (KamLAND Collaboration), "First Results from KamLAND: Evidence for Reactor Antineutrino Disap- 
pearance," [arXiv:hep-ex/0212021]. 

[8] M. C. Gonzalez-Garcia, "Theory of Neutrino Masses and Mixing," Plenary report at the 2003 International Conference on 
High Energy Physics, Amsterdam, [arXiv:hep-ph/0210359]. 

[9] M. Apollonio, et al. (CHOOZ Collaboration), [arXiv:hep-ex/0301017]. 
[10] F. Boehm, et al., Phys. Rev. D64, 112001 (2001) [arXiv:hep-ex/0107009]. 

[11] J. G. Learned and S. Pakvasa, Astropart. Phys. 3, 267 (1995); H. Athar, M. Jezabek, and O. Yasuda, Phys. Rev. D62, 
103007 (2000). Osamu Yasuda, "Neutrino oscillations in high energy cosmic neutrino flux," [arXiv:hep-ph/0005135]; G. J. 
Gounaris and G. Moultaka, "The Flavor Distribution of Cosmic Neutrinos," [arXiv:hep-ph/0212110]. 

[12] J. P. Rachen and P. Meszaros, Phys. Rev. D58, 123005 (1998). 

[13] V. D. Barger, D. Marfatia, and B. P. Wood, Phys. Lett. B498, 53 (2001); see also http://www.awa.tohoku.ac.jp/KamLAND/ 
[14] http://www-numi.fnal.gov 



9 



[15] http:/ /proj-cngs.web.cern.ch/proj-cngs/ 

[16] G. Barenboim, L. Borissov, J. Lykken, and A. Yu. Smirnov, J HEP 0210, 001 (2002) [arXiv:hep-ph/0108199]; G. Barenboim, 

L. Borissov, and J. Lykken, Phys. Lett. B534, 106 (2002) [arXiv:hep-ph/0201080]. 
[17] G. Barenboim, L. Borissov, and J. Lykken, "CPT Violating Neutrinos in the Light of KamLAND," [arXiv:hep-ph/0212116]. 
[18] J. F. Beacom, N. F. Bell, D. Hooper, S. Pakvasa, and T. J. Weiler, "Decay of High-Energy Astrophysical Neutrinos," 

[arXiv:hep-ph/0211305]. 
[19] J. F. Beacom and N. F. Bell, Phys. Rev. D65, 113009 (2002) [arXiv:hep-ph/0204111]. 

[20] The interaction length of a 1-TeV neutrino is about 2.5 million kilometers of water, or 250 kton/cm 2 , whereas high-energy 

photons are blocked by a few hundred g/cm 2 . 
[21] We use the symbols $ and <p to denote the sum of neutrinos and antineutrinos except when we discuss implications of 

possible CPT violation, where $, refer to neutrinos, and $, (p to antineutrinos. 
[22] Muon cooling within the source, as discussed in Ref. |T3l . would deplete the source flux of We have no candidate 

mechanism to produce an ultrahigh-energy flux like $?. In view of the current puzzlement over the origins of ultrahigh- 

energy cosmic rays, we think it prudent to keep an open mind about the range of possibilities. 
[23] The barred quantities refer to antineutrinos. 
[24] A similar phenomenon is familiar for cosmic-ray muons. 
[25] Our assumption that Tz/mz = T2/7712 = r/m is also a special case. 



